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Topological superconductors are gapped superconductors with gapless and topologically robust 
quasiparticles propagating on the boundary. In this paper, we present a topological field theory 
description of three-dimensional time-reversal invariant topological superconductors. In our theory 
the topological superconductor is characterized by a topological coupling between the electromag- 
netic field and the superconducting phase fluctuation, which has the same form as the coupling of 
"axions" with an Abelian gauge fleld. As a physical consequence of our theory, we predict the level 
crossing induced by the crossing of special "chiral" vortex lines, which can be realized by considering 
s-wave superconductors in proximity with the topological superconductor. Our theory can also be 
generalized to the coupling with a gravitational fleld. 



00 



o 
o 

I 



I 

O 

o 



> 
O 

o 



X 



I. INTRODUCTION 

Since the discovery of the quantum Hall effect in the 
1980'sii^, the study of topological states of matter (TSM) 
has been an active field of research in modern condensed 
matter physics. TSM arc those states of matter which 
are distinguished from each other by some topological 
properties rather than more conventional properties such 
as symmetries preserved or broken by the order param- 
eter. Recently, wide classes of TSM's known as topo- 
logical insulators (TI) and topological superconductors 
(TSC) have been proposed, and TTs have been experi- 
mentally realized in various materials^"—. TPs are insu- 
lators with a bulk gap and gapless edge/surface states 
propagating on the boundary which are robust against 
perturbations within the given symmetry class. Sim- 
ilarly, TSC's are fully gapped superconductors with a 
bulk gap for the quasi-particle excitations, and topologi- 
cally protected gapless quasi-particle states propagating 
on the boundary. The first example of topological super- 
conductor is the (p + ip)-wave superconductor of spinless 
fermions in two-dimensions^, which has a chiral Majo- 
rana fcrmion state propagating on the edge. (The Majo- 
rana fermion of the edge state refers to the fact that the 
quasiparticle is its own anti-particle.) More recently, a 
new class of TSC was proposed in three-dimensional (3d) 
time-reversal invariant (TRI) superconductors^i"— , which 
has spin-triplet pairing in the bulk, and two-dimensional 
masslcss Majorana surface states with linear "relativis- 
tic" dispersion. The "^Hc B phase was proposed as an 
example of the 3d TRI TSC (more rigorously, a topolog- 
ical superfluid)iii. Solid state candidates of 3d TSC have 
also been propose d^^i^^ . A recent tunneling experimeniii^ 
supports the theoretical proposal of Cu doped Bi2Se3 as 
a TSCi^, although more experiments are needed to reach 
a conclusion. 

Without interactions, TI and TSC can be described by 
topological band theories^, which characterize a class of 
TI or TSC by a topological invariant of the band Hamil- 
tonian. The classification of non-interacting TI and TSC 
is understood systematically in generic dimensions^i d^i^^ . 
However, in reality all electron systems are interacting. 



so that the physically interesting TSM are those which 
are robust even with interaction. In general it is difficult 
to directly study the topological classification of gapped 
interacting Hamiltonians or their ground states, except in 
some special cases such as in one dimcnsioni^"— , and in 
some special classes of models in higher dimensions^Sri^l. 
A general approach to characterize interacting TSM's is 
by writing down the possible topological response theo- 
ries, which describe some observable physical properties 
of the system, and contain a "topological order parame- 
ter" that is required to be quantized by general princi- 
ples. For example, the quantum Hall effect can be charac- 
terized by Chern-Simons theories^^"— . The topological 
order parameter is the coeflScient of the Chern-Simons 
theory which has the physical meaning of Hall conduc- 
tivity and is quantized due to charge quantization. In 
the case of 3d TRI TI, the topological response theory 
is described by a term in the electromagnetic responsei^ 
/ d'^xF^i/Fr^r^^""'^ , in which 6 is quantized to be 
or TT mod 2-k if the system is time-reversal invariant. This 
is an example of a topological response theory in which 
the quantization of the topological order parameter de- 
pends on discrete symmetry (time-reversal symmetry in 
this case). The physical consequence of this topologi- 
cal response theory is the topological magneto-electric 
effec t^^i^^ . Such a topological term has the same form as 
the coupling of axions with gauge field proposed in high- 
energy physics^Si^. However, in a TI, 6' is a constant 
determined by the bulk topology rather than a dynami- 
cal field. 

It is therefore natural to ask what is the topologi- 
cal response theory description of TSC. In a supercon- 
ductor the electromagnetic field is screened, and the 
edge/surface states of TSC are quasi-particles without 
a well-defined charge. Therefore the electromagnetic re- 
sponse has been considered as improper for defining the 
topological response theory. For this reason, the topo- 
logical BF theory has been proposed to describe the 
chiral topological superconductor—. Alternatively, the 
coupling to gravitational field and quantized thermal re- 
sponse have been proposed^'22'^. Similar to the axionic 
coupling to the electromagnetic field in 3d TI, a term 
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proportional to the Pontryagin invariant of the Ricmann 
curvature is proposed to characterize the gravitational 
response of 3d TSC, with coefficient 9 = or n mod 
27p32j33 However, Ref. [Hj pointed out that the coef- 
ficient 9 of the Pontryagin invariant is a U{1) phase, so 
that under the condition of time-reversal symmetry there 
are only two physically different values of 9, leading to 
a Z2 classification of the TSC rather than the integer 
classification in the non-interacting case. Ref. [s^] pro- 
posed an alternative approach of defining a topological 
response theory by considering proximity effects of TSC 
to trivial s-wave superconductors. On the surface of a 
TRI TSC, two thin films of trivial s-wave superconduc- 
tors are coupled to the TSC by the proximity effect and 
also coupled with each other through Josephson coupling. 
When the two films of the s-wave superconductors have a 
TT phase difference, the number of chiral modes propagat- 
ing along the Josephson junction provides a definition of 
the integer valued topological invariant of the bulk TSC. 
This proposal provides a physically operational definition 
of the topological invariant, which answers the question 
of whether the integer valued topological invariant of 3d 
TRI TSC is robust upon electron interaction. However, 
the s-wave proximity effect is proposed as an experimen- 
tal setting and is not formulated in terms of a topological 
field theory description. 

In this paper, we develop the topological field theory 
description of 3d TSC by considering the TSC coupled to 
both the electromagnetic field and the superconducting 
phase fluctuations. We demonstrate that the electromag- 
netic response of the TSC contains a topological term 

^ _ -(T^ Cii9i , . 

i 

with i labeling all Fermi surfaces and 9i the supercon- 
ducting phase of each Fermi surface. Cu denotes the 
Chern number describing the spin helicity on each Fermi 
surface. We dub this theory the "axion field theory" of 
TSC since it has the form of the coupling between axion 
and Abelian gauge theory^^i^. The details of the deriva- 
tion and physical consequences will be presented in the 
later part of this paper. The form of this topological field 
theory appears similar to that of TI, but there are two 
essential differences: i) The phases 9i are dynamical vari- 
ables describing the superconducting phases, which are 
thus coupled to the gauge field by gauge coupling, 
while in TI the phase is a constant topological order pa- 
rameter; ii) In a superconductor there is a Higgs term 
for the gauge field A^, so that the gauge field is mas- 
sive, while in TI the gauge field remains gapless in the 
long wavelength limit. Due to these two important differ- 
ences, the physical consequences of this topological term 
are very different from those of the TI case. Due to the 
Higgs term, the electromagnetic field vanishes in the su- 
perconductor except in vortex cores (if the superconduc- 
tor is type II) . For a superconductor with multiple Fermi 
surfaces, we can consider vortices of only one Fermi sur- 
face, which we call chiral vortices since each Fermi surface 



with a nontrivial Chern number is equivalent to a Weyl 
fermion with chirality. We show that the topological term 
([1]) describes an anomaly when chiral vortices cross each 
other, which can be interpreted as the change of ground 
state fermion number parity caused by level crossing be- 
tween quasi-particle energy levels. This is an analog of 
the Z2 Wittcn anomaly in a system of odd flavors of Weyl 
fermions coupled to an SU{2) gaugcfield^i. The s-wave 
proximity effect proposed in Ref. [3^ is shown to be a 
physical realization of the chiral vortices. Therefore our 
topological field theory correctly describe the proposed 
s-wave proximity effect^. 

The rest of the paper is organized as follows. The sin- 
gle particle description of 3d TRI TSC is reviewed in 
Sec. ini which provides a starting point of our discussion. 
In Sec. mil we present the derivation of the axion field 
theory by introducing a mapping from the 3d TSC to a 
four-dimensional TSC. In Sec. II VI we discuss the physical 
consequences of the axion field theory, including the chi- 
ral vortex linking effect and the s-wave proximity effect. 
Finally, Sec. |V] is devoted to the summary and further 
discussions. 



II. THE (3 + 1)-D TOPOLOGICAL 
SUPERCONDUCTORS 

A. The topological invariant of (3 + l)-d topological 
superconductors 

As a starting point of our theory, we review the single 
particle description of the (3-1- l)-d topological supercon- 
ductors. At the level of Bogoliubov-de Gennes (BdG) 
theory, the topological superconductors are described by 
the generic BdG Hamiltonian^. 




with Ck a A^-component vector for a system with A^ 
bands, and /ik, Ak N x N matrices. Time reversal sym- 
metry acts as 

r(ck) = rc_k, T(cLk) =r*4 (3) 

with T the time-reversal matrix satisfying = T* — 
— T, = —1. The time- reversal invariance requirement 
to the Hamiltonian is 

T^hiT = /i-k, T^A^r = -A_k (4) 

By defining Qk = + *Ak, one can prove that Qk is 
non-singular for gapped superconductors. The integer- 
valued topological invariant describing three-dimensional 
topological superconductor is the winding number of the 
map from the Brillouin zone torus to the special linear 
group SL{N, C) defined by Qk- 

For the purpose of this paper, it is helpful to start 
from an alternative formula of the topological invariant 
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proposed in Ref. [35[ in the weak pairing limit. When the 
pairing is weak compared to the kinetic energy /i^, only 
the pairing near Fermi surfaces of the kinetic energy /ik 
is important for the topological properties of the system. 
For each Fermi surface I]„ , the momentum dependence of 
the states jyik) on the fcrmi surface defines a Berry con- 
nection fli = —i {nk.\ dki |fik). which has a Chcrn number 
defined by 

Cm^^j dn'^d^a, - d,a,) (5) 

If one defines the pairing order parameter near the Fermi 
surface by 

A„k = r((nk|)Ak|7ik) (6) 

it can be proved that A„k is real due to time-reversal 
symmetry. The choice of r((nk|) as the wavefunction at 
momentum — k avoids the ambiguity in the gauge choice 
of the wavefunction, and makes the statement of a "real" 
pairing A„k meaningful, which is only possible because of 
the time-reversal symmetry. The gap of the BdG Hamil- 
tonian near the Fermi surface is |A„k|, so that in a fully 
gapped superconductor, A„k cannot change sign as a 
function of k on each Fermi surface. Consequently, in 
a fully gapped TRI superconductor, the sign of the pair- 
ing is definite in each Fermi surface. It was proved in 
Ref. [ssj that the topological invariant is given by the 
following sum of Fermi surface Chcrn numbers weighted 
by the sign of the pairing: 

^^Ci„sgn(A„k) (7) 

n 

As an example, we can consider the following two-band 
model defined by 

k2 

/ik ~ /i + aa ■ k 

2m 

Ak = iAoCTj^cr • k (8) 

in Eq. The time-reversal transformation is defined 
by 7" = icTy. Up to a basis transformation, the Hamilto- 
nian for a = describes the Hc'^ BW phase^Ti^i^. We add 
a term a to lift the degeneracy between the two Fermi 
surfaces, which is convenient for defining the pairing or- 
der parameter on the Fermi surface. We would like to 
note that this splitting is just for convenience and is not 
physically required. 

The two Fermi surfaces are spherical with the Fermi 
momentum 

kp± — =pma -|- •\/ m^a^ + 2mfj.. (9) 

The states around the two Fermi surfaces are defined by 
the two eigenstates of ct • k, denoted by cr • k |k±) = 
±|k| |k±). It can be verified that the two Fermi surfaces 
with the opposite spin helicities have opposite Chern 
numbers and also opposite signs of the pairing order pa- 
rameter, so that the topological invariant is given by 

iV = i [1 - (-1)] sgnAo = sgnAo (10) 



B. A global ambiguity in the definition of the 
topological invariant 

Before proceeding to the next section, we would like to 
make a comment about the definition of the topological 
invariant ([7]). The pairing order parameter A„k is defined 
in Eq. (j6|), which depends on the definition of the time- 
reversal transformation in Eq. ([3]). If we substitute T 
by — T in Eq. all properties of T are still preserved, 
and the time-reversal condition for Hamiltonians given 
by Eq. ^ also remains valid. Therefore a time-reversal 
invariant Hamiltonian remains so in the new definition of 
time-reversal. However, A„k defined in Eq. ^ changes 
sign for all bands. Therefore the topological invariant 
N is transformed to —N. Since changing the sign of T 
is equivalent to following the time-reversal transforma- 
tion by a global phase rotation Ck — > — Ck, there is no 
physical way to distinguish the two definitions of time- 
reversal symmetry. Consequently, topological supercon- 
ductors with topological invariant and —N arc physi- 
cally equivalent. However, such an ambiguity only occurs 
for a global sign of the topological invariant and does not 
affect any physical consequence of the topological invari- 
ant such as surface states at the interface between two 
superconductors with different topological invariant. It 
also docs not affect the topological field theory and topo- 
logical defects discussed in the rest of this paper. 

III. (4 + l)-DIMENSIONAL REGULARIZATION 
AND THE DERIVATION OF THE AXION FIELD 
THEORY 

With the background prepared in the last section, now 
we show how the theory reviewed above is related to an- 
other topological state of matter, the (4-1- l)-dimensional 
time-reversal invariant topological insulatorsi^i^. Such 
a relation will be essential in the axion field theory that 
we will derive. 

A. (4 + l)-dimensional model of the (3 + l)-d TSC 

From the Fcrmi surface formula, it can be seen that 
Fermi surfaces with nontrivial Chern number are essen- 
tial for topological superconductivity. On the other hand, 
a Fermi surface with non-trivial Chern number Ci = 1 
is topologically equivalent to a Weyl fermion (with the 
fermi level generally away from the Dirac point). For ex- 
ample, in the minimal model the low energy effective 
theory near each Fermi surface is a Weyl fermion since it 
has a fixed spin hclicity a ■ k/|fc|. The Fermi surface at 
kp^ has spin hclicity +1 with a uniform pairing Aokp^. 
Thus the low energy states around the spin helicity -1-1 
Fcrmi surface in the weak pairing limit is equivalent to 
the following Weyl fermion theory: 

Hr = ^ V'flk ["F+cr ■ k - 

k 
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(11) 



with = VF+kF+ and vpj^ is the Fermi velocity at this 
Fermi surface, h.c. stands for Hermitian conjugate of the 
pairing term. It should be noticed that a uniform pair- 
ing is used here to reproduce the pairing near the Fermi 
surface in the model (|8]), which is possible because the 
structure of the pairing is determined by time-reversal 
symmetry, such that the topologically trivial and non- 
trivial pairings have the same form when restricted to 
one Fermi surface. The main difference between topolog- 
ically trivial and nontrivial pairings is the relative sign 
between the pairing order parameters at the two Fermi 
surfaces. The other Fermi surface with spin helicity —1 
can be represented similarly by a Weyl Fermion with op- 
posite helicity: 



H, 



Lk 



-\\Y.^okF+i^l-\.^^y^l\. + h.c^ (12) 

The Hamiltonian Hl + Hr correctly reproduces the 
low energy behavior of the minimal model ([5]) in the 
weak pairing limit. More generically, one can represent 
each Fermi surface with Chern number 1 (—1) by a Weyl 
fermion with right-hand (left-hand) helicity. The Fermi 
surfaces with higher Chern number can be similarly rep- 
resented by multiple channels of Weyl fermions. 

The advantage of the Weyl Fermion representation is 
its relation to the (3-|-l)-d topological superconductor to 
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(4-1- l)-d topological insulators. As is proved in Ref. 
time-reversal invariant (4 -f l)-d insulators are classified 
by an integer valued topological invariant, the second 
Chern number 



Co = 



1 
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(13) 



with /i. 



djUi 



+ i \ai,aj 



—i (nk| dki |mk) the U{N) Berry phase gauge field for a 
system with N bands occupied. |nk) labels the occupied 
bands. Coupled with a charge U{1) gauge field A^, the 
system with Chern number C2 has a topological response 
described by the (4 + l)-d Chern-Simons term 



^CS [Aa] 



C2 

247r2 



(14) 



with a, 6,.. denoting the (4 + l)-d space-time indices 
0,1,2,3,4. In a geometry with open boundary, the 
(4 -I- l)-dimensional topological insulators have surface 
states on the (3 -|- l)-d boundary described by Weyl 
fermions. The net helicity of Weyl fermions is given 



by 



iV_ 



C2, with N± the number of right- 



handed/left- handed Weyl fermions. The charge trans- 
port determined by the Chern-Simons term (jl4p 
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^abcde tt^ r? 
e ^bc't'de 



(15) 





FIG. 1: Representation of the 3-|-ld topological superconduc- 
tor on the spatial manifold M3 as the boundary state of a 
4+ld topological insulator defined on M3 x /. 



correctly reproduces the axial anomaly of the Weyl 
fermions on the boundary^. 

Now consider a (4 -I- l)-d topological insulator with 
Chern number C2 = 1 on a slab geometry x I (with 
the torus and / the interval), as shown in Fig. [T](a). 
For example a representative Hamiltonian in this class is 
the lattice model introduced in Ref. [l4l: 
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sin hr 



COS ki 



Ck 



(16) 



with F", a = 0, 1, .., 4 Hermitian matrices satisfying the 
Clifford algebra {r°,r''} = 26"''. For -4 < m < -2 
the system is a topological insulator with Chern num- 
ber C2 = 1, and the surface state on each boundary of 
X I consists of a single copy of Weyl fermion, with 
opposite helicity on the two boundaries. The system is 
time- reversal invariant, with the time-reversal transfor- 
mation matrix T satisfying T~^r*T = — F^, i = 1, 2, 3, 4, 

We denote the periodic dimensions as xi^2.3 S 
[0, ^1,2.3] and the dimension with open boundary con- 
ditions as Xi G [0, L4]. In the limit of finite L4 and 
-^1,2,3 00, the (4 -I- l)-d system can be considered as a 
(3 + l)-d system with periodic boundary conditions, and 
a large number (4L4) of bands per site. However, the low 
energy effective theory of this system only consists of two 
massless Weyl fermions with opposite helicity. 

Now we relate this model to the (3 -I- l)-d topological 
superconductor by considering a superconducting pair- 
ing term at each boundary of a(4 + l)-d system. Phys- 
ically, one can consider the geometry shown in Fig. ([1]) 
(b), with the topological insulator sandwiched between 
two thin layers of (3 + l)-dimensional ordinary s-wave 
superconductors. Denoting the phases of the two super- 
conductors hy 9l and 9i^, the low energy effective theory 
of this sandwich structure is given by 

H ^ [i'R)^a- ■ kV'flk - ^IkCr ■ k^pLk 

k 
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iayip 



t 



h.C. 



(17) 



Compare this effective theory with Eq. ((TT|) and (IT2]) . 
one can see clearly that the SC-TI-SC sandwich struc- 
ture is topologically equivalent to a (3 + l)-d topologi- 
cal superconductor described by the Hamiltonian (|5]) if 
Ol = 0, Or = t:. 

In a similar way, the topological superconductor with 
higher winding number TV can also be represented by a 
similar (4 + l)-d "sandwich" model with the bulk Chern 
number C2 ~ N, and two surface superconductors with 
opposite sign of pairing. The surface states on each 
boundary consist of N Weyl fermions, each of which has 
a Fermi surface with first Chern number Ci = 1. If the 
two surfaces have opposite sign of pairing, one can verify 
that the topological invariant of the 3d TSC is N. 



B. Axion field theory 

The mapping of the (3 + l)-d TSC to the (4 + 1)- 
d model discussed above may seem uneconomic at first 
sight, since the (4-t- l)-d model contains many more high 
energy degrees of freedom. However, the (4-1- l)-d model 
has the advantage of spatially separating the left- and 
right-handed Weyl fermions to different surfaces, which 
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is convenient for deriving the effective field theory de- 
scription of the TSC. To obtain the effective field theory, 
we consider a uniform superconducting pairing term on 
the boundaries 



Hsc = ^Aoe-^«^("-)c^(x^,a;4 = 0)rc(x^,a;4 = 0) + /i.c. 

+ Aoe"''«("'')c^(a:^, a;4 = L4)rc{x^, X4 - L4) 
+h.c. (18) 



Here the two terms are pairing of electrons on the left and 
right boundary defined by the 4-th coordinate X4 = and 
L4 respectively. We consider the pairing amplitude Aq as 
a constant, and keep the phase fluctuation 9r(^^(^x^) as 
generic fields depending on the boundary (3-|-l)-d space- 
time coordinates x^, ^ = 0,1, 2, 3. The physical electro- 
magnetic field in the (3 -I- l)-d space-time is represented 
by a U{1) gauge field defined in the (4-1- l)-d space-time, 
with the restriction that ^4 = and = A^j^^x^) only 
depends on the (3 + l)-d coordinates x^. The fermions 
are minimally coupled to the gauge field. The effective 
action of the boson fields 6b,0l,A^ is defined by inte- 
grating out the fermion in the action: 



exp(iS'cff [A^,6'L,6'i?]) = J DcDcexp {iS [c,c, Afj,,9L, Or]) 
S[c,c,A^,9l,0r] = J 



Cj {dt - iAo) c, - i?buik [c, c, Ai] - iJsc [c, c, Ol,Or] 

. i 

I 



(19) 



with -ffbuik [c, c, Ai] defined by the Hamiltonian of topo- 
logical insulator given in Eq. (|16|) , and Hsc the boundary 
pairing terms P^ . 

To compute the effective action it is helpful to notice 
the gauge invariance of the action S [c, c, A^, 0^, Or] un- 
der the following gauge transformation: 

c{Xa) ^ ciXaV^^^^ 
Aa{Xb) Aa{xb)+daip 

Or{x^j,) Or{x^) + 2(p{x^,X4 = 0) 

OUx^) ^ 9Lix^) + 2ip{x^,Xi^ U) (20) 

Here Xa,a = 0,1, 2, 3, 4 denotes the (4 + l)-d coordinates. 
In particular if wc make the special choice 



1 

2i^ 



[0LiL4 - X4) + ORX4] 



(21) 



Ol , Or are canceled by the gauge transformation and the 
gauge field is transformed to Aa — Aa + da^p- This 



way, the dependence of the effective action on the pair- 
ing phases 0l,9r has been absorbed into the dependence 
to the C/(l) gauge field, and the U{1) gauge field now 
obtains the 4-th component A^ ~ d^Lp = {Ol — Or) /2L4. 

Since the pairing terms are only introduced on the 
boundary, the bulk remains a (4 -I- l)-d topological in- 
sulator coupled with the gauge field Aa- Consequently, 
the effective action of the gauge field still contains the 
the Chern-Simons term given in Eq. (|14p with coefficient 
C2 = 1. Written in the physical degrees of freedom 6'i?,(L) 
and A^, we have 



S'cff = Scs 
1 



1 



247r2 



d^xe-'-^'^AadbA^d^A, 



32^2 / d'xe^^'^'^^^-^^^F^^F^r. (22) 

Thus we see that the Chern-Simons term on the slab 
geometry is reduced to the (3 -|- l)-d axion term with 
the U{1) field di — Or coupled to the gauge field in an 
axionic coupling. Besides the topological term (|22p . we 
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also have the ordinary terms in a superconductor, in- 
cluding Higgs term and the higher order terms such as 
the Maxwell term. Since the two boundaries represent 
the two Fermi surfaces in the physical superconductor, 
the two superconducting phases 0^, 6*^ are generically 



not decoupled since there is only a global U{1) symme- 
try. Thus Josephson type coupling terms cos(6'i — are 
also allowed. With these considerations, we finally obtain 
the following effective action of the (3 + l)-d topological 
superconductor : 



5'off 



64^2 ^M-^- 



462 '^" 



2A^)'-f Jcos {0l-Or) 



(23) 



The effective action ((23)) is the central result of this 
work. It describes the response of the fermions to super- 
conducting phase fluctuations Oj^^Oj^ of the two Fermi 
surfaces, and to the electromagnetic field. It should be 
noticed that the effective action applies to both topologi- 
cal and trivial superconductors depending on the ground 
state values of OljOr- As has been discussed in Sec. II, 
= 0, 0R = TT corresponds to a topological supercon- 
ductor, and Ol = Or = corresponds to a trivial super- 
conductor. In the effective action (|23|) . the ground state 
value oi 9l — 9r can be controlled by the sign of J. J > 
(J < 0) describes topological (trivial) superconductor. 
Since this effective theory is obtained by integrating out 
fermions, it only applies to the energy scale below the 



BCS gap E < \Al\, |Ai?|. Therefore the coupling be- 
tween the two Fermi surfaces must be weak enough such 
that the mass of the mode 9l — Sr m? = 2J / p <^ IAl^^P, 
in order for the effective theory to be meaningful. 

The approach of (4 -|- l)-d regularization can be gen- 
eralized to generic TSC with multiple Fermi surfaces. In 
the weak pairing limit, a superconducting phase 9i can 
be defined on each Fermi surface. If a Fermi surface has 
the Chern number Cii, it is equivalent to Cu copies of 
the Weyl fermion studied above. The axion term in the 
electromagnetic response is given by a sum over all Fermi 
surfaces. Thus we obtain the following generic effective 
action: 



5'off = / d^:i 



1 

64^^' 



(24) 



It is worthwhile to note that a time-reversal transforma- 
tion acts by 9i -9i, Aq ^ Aq, A^^^y^z -Ax,y,z, and 
the effective action is invariant, as expected. The Joseph- 
son coupling between the superconducting phases of dif- 
ferent Fermi surfaces is always in the form of cos{6i — 6j) 
as is required by time-reversal symmetry. The ground 
state value of 9i (up to an overall constant) is deter- 
mined by the Josephson couplings Jij. In general, it is 
possible that the coupling is frustrated and the ground 
state spontaneously breaks time-reversal symmetry, al- 
though we are mainly interested in the situation with 
time-reversal symmetry preserved, with 9i = or tt for 
each i . In that case the topological invariant ([7]) is given 
by 

i 

To summarize this section, we have obtained the 
generic effective action of (3 -|- l)-d superconduc- 



tors coupled with the electromagnetic field and supercon- 
ducting phase fluctuations. The most important term in 
the action is the axionic coupling between the supercon- 
ducting phase fluctuations and the electromagnetic field, 
which is a consequence of the nontrivial Chern number 
of the Fermi surface. The effective action describes both 
topological trivial and nontrivial superconductors. Al- 
though not explicit in the effective action, the topological 
invariant N is determined by the Chern numbers Cu and 
the ground state value of Oi, which is thus determined for 
a given effective action (|24l) . 



IV. TOPOLOGICAL EFFECTS DESCRIBED BY 
THE AXION FIELD THEORY 

One of the main advantages of an effective field the- 
ory description is that it directly describes observable 
topological response properties of the topological state of 
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matter. Having obtained the effective field theory (|241) . 
we now discuss the topological effects described by this 
theory. 



A. Anomaly and chiral vortex lines 



Naively, the effective action (|24p is dominated by the 
Higgs terms — 2Af^) , which leads to Meissner effect 
and screens the electromagnetic field in the superconduc- 
tor. With zero field strength F^^, in the superconductor, 
the topological term will have no physical consequence. 
However, the situation becomes nontrivial when vortex 
lines are considered in a type-H superconductor. Around 
an ordinary Abrikosov vortex line, the phases Oi on all 
Fermi surfaces have the same winding number, so that 
the phase combination CuOi in the topological term 
has no winding. To see the nontrivial consequence of the 
topological term, it is essential to consider a "chiral vor- 
tex line" where only some of the fl^'s have a vorticity. In 
the following we will study the anomaly at presence of 
such chiral vortex lines, as a consequence of the topolog- 
ical term. 

We start by considering the equation of motion of the 
gauge field determined by the action (p4)) : 

= ^ (26) 



The charge current is defined as 



= 2^p, {d^e,-2A, 



(27) 



which should be conserved due to the global U{1) sym- 
metry. However, with the F A F term in the action ([24| . 
there is an anomaly in charge conservation. According 
to Eqs. (1211) and (I^, we obtain 



i 



with 



In 



(28) 



(29) 



the antisymmetric tensor field for the vortex current. For 
example, if there is only a vortex line for 6i along the 
z direction at coordinate a; = j/ = 0, with the Chern 
number Cn = 1, we have 



(30) 



with other components vanishing. The anomaly is given 
by 



J J vortex 



(31) 



which is half of the Chern number of the electric field 
F^jj in the zt plane. Physically, such an anomaly is re- 
lated to the fact that there is a 1 + 1-D Majorana-Weyl 
fermion propagating along the vortex line^. The anoma- 
lous term above obtained from the effective theory can- 
cels the anomaly of the Majorana-Weyl fermion. By com- 
parison, in the more familiar case of a 1-1- Id Weyl fermion 
coupled to gauge field A^, the anomaly is given by 



47r 



(32) 



with a,b ~ 0,1. For an electric field with fiux 27r, 
J cPxd°'ja = 1 which means that one additional Weyl 
fermion appears in the 1 -I- 1-d system. Physically, one 
can view the Weyl fermion as an edge state of the quan- 
tum Hall system, and the anomaly is interpreted as a 
charge pumping from the bulk to the boundary due to 
the Hall currenti^ Since a Majorana-Weyl fermion car- 
ries half the degree of freedom of a Weyl fermion, the 
anomaly ((3T|) is also given by half of that of the Weyl 
fermion. 

To obtain a more physical understanding of the 
anomaly equation, the motion of multiple vortex lines 
needs to be considered, since the electricmagnetic field 
inducing the anomaly is only present in the superconduc- 
tor at vortex lines. For simplicity, consider the = 1 
topological superconductor with two Fermi surfaces sat- 
isfying 6l — 6ii = TTas an example. According to the 
Higgs term, away from the vortex lines we have 

PL{d^0L-2A^)+PRid^eR-2A^) = (33) 

so that around a chiral vortex line of 6^, the fiux is 



A-dl = 



1 



PL 



PL + PR 



pr) 

-TT 



L + Pfl,V6ifl) • d\ 

(34) 



in which the loop integral is taken around a loop enclos- 
ing the chiral vortex line. Therefore the chiral vortex 



line carries a fiux of 



PL he 
Pl+Pr 4e 



in Gauss units. Similarly, 



the chiral vortex line of Or carries a flux of — — 
Consider the configuration illustrated in Fig. 



l+Pr 4e ■ 

H with 4 



chiral vortex lines, 2 for each of 9^ and Or, along the y 
direction, and 2 ordinary vortex lines along the z direc- 
tion. Periodic boundary conditions are assumed along all 
three directions x, y, z. It should be noted that in such a 
compact manifold, the fiux of the gauge field in each 
plane xy, yz, zx must be quantized in unit of 2tt, which is 
why the number of chiral vortex lines must be a multiple 
of 4, and that of ordinary vortex lines must be even. 

Now consider the motion of z vortex lines across the 
chiral ones, as is illustrated in Fig. O When both of the 
z vortex lines move from the left to the right of the chiral 
ones, due to periodic boundary conditions, they can re- 
turn to their original position, so that such a process is a 
periodic evolution of the system. In the following we will 
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show that such a vortex motion corresponds to an instan- 
ton process of i^^^ with nontrivial / d'^xe^'^'^'^ F^^F^t- 

Since the field strength is only nonzero in the vortex 
cores, the term e^^'^'^ F^^F^r vanishes except when two 
vortex lines cross each other. The field strength corre- 
sponds to a static y chiral vortex line at position (xo, zo) 
is 

F^l^ = \5{x - xo)5{z - zo) (35) 

with other components vanishing. Similarly for an ordi- 
nary vortex line along the z direction with coordinates 
{xQ^t) = vt^yo), the field strength is 

i^i^) = -KSix ~ vt)5{y ~ y^) 

F,f = -vnSix - vt)Siy - yo) (36) 

Therefore the contribution of the two vortex lines to the 
topological term is 

= Att'^v5{x - vt)S{y - yo)S{x - xo)S{z - zq) 
= 47r^(5(t - —)S{y - yo)5{x - xo)S{z - zq) 

V 

= 4tt^6\x^ - xo^) (37) 

with (5^ denoting the 4-dimensional 5 function and xof^ = 
{xq/v, xo,yo, Zq) the crossing point in space-time. In the 
process shown in Fig. [51 there are 8 such vortex crossing 
points. If we denote the location of the 8 crossing points 
by x^j, n = 1,2,..., 8, the topological term for such a 
process is 

Jtopo — U X ^tJ.iy-rcrT 
1 ^ 

= Ye^^^^^'^"^"^^^"'"^^ ^^^^ 

n— 1 

To understand the implication of such a term, consider 
the transformation 9l (x) — >■ 9l (x) + 27r with Oji invari- 
ant. The change of the topological term is AS'topo — tt, 
so that the partition function changes sign Z — >■ —Z. 
This seems to be contradictory with the fact that 0l is 
the phase of the order parameter and is periodic with the 
period of 2tt. However, this is exactly the consequence of 
the anomaly in the Majorana fermion system. The phase 
rotation 9l 9 + 2tt corresponds to a phase rotation of tt 
for electrons iPl{x) — > In such a transformation, 

a many-body state obtains a phase (—1)^^ with Nl the 
number of left-handed fermions. (—1)^^ is the fermion 
number parity of the left-handed fermions. In the adia- 
batic limit, the propagator e''^'°p° can be considered as a 
propagator of the ground state back to itself: 

g^Stopo ^ i^Q\ ^-'S-Z \G) (39) 

The fact that e*'^'°p° — ~e^^^°^° in the global transfor- 
mation —i^L can be interpreted as the initial and 



final states having opposite fermion number parity for 
the left-handed fermion. Similarly, the initial and final 
state also have opposite fermion number parity for the 
right-handed fermion. This effect can be understood as 
a residual effect of the anomaly in Weyl fermions. If we 
consider the normal state of the topological supercon- 
ductor, it is topologically equivalent to a pair of Wcyl 
fermions with opposite chirality. In a background field 
F^i^, the Weyl fermions have the anomaly 

9,31 = -« = ^^"''^'^M.^- (40) 

In an instanton configuration with C2 = 
32^e'^''--F^,F,, = 1, we have ! d^xd^j^ = 
— J d^xd^j^^ = 1, which means one left-handed Weyl 
fermion becomes right-handed after the instanton pro- 
cess. When we consider superconductivity, the fermion 
number conservation is broken, but the fermion number 
parity conservation is still preserved. Consequently the 
anomaly is still well-defined for instanton processes with 
odd C2. The vortex line crossing process we discussed 
above is an example of a configuration with C2 = 1. 

More microscopically, such an anomaly can be related 
to some spectral fiow of the Majorana fermions, as is il- 
lustrated in Fig. [21 This is similar to the spectral flow 
picture of the Z2 Wittcn anomaly^!. Each chiral vor- 
tex line hosts a (1 -(- l)-d Majorana- Weyl fermion, which 
can be seen by explicitly solving the Majorana fermion 
spectrum with a vortex line of the mass term. This is 
the Majorana version of the axion string. When peri- 
odic boundary conditions are considered for the chiral 
vortex lines along the y direction, whether or not each 
chiral vortex line has a Majorana zero mode at exactly 
zero energy depends on the boundary condition of the 
Majorana- Weyl fermion. We start from the configuration 
in Fig. [2I (a) and assume that the horizontal vortex lines 
have anti-periodic boundary conditions and thus have no 
zero modes. When the perpendicular vortex line crosses 
the first horizontal one, the boundary condition for the 
Majorana- Weyl fermion along the horizontal chiral vor- 
tex line is changed, so that a pair of Majorana zero modes 
appear on the two vortex lines which have just crossed 
each other. Subsequently, when the perpendicular vortex 
line crosses the second chiral vortex line, there are two 
Majorana zero modes on the two horizontal vortex lines. 
When the second perpendicular vortex line crosses them, 
the boundary condition is changed back to anti-periodic 
and thus the zero modes disappear. The spectral flow in 
the whole process is summarized in Fig. [21(e) which con- 
tains a level crossing between electron and hole states. 
Due to this level crossing, one would conclude that the 
ground state on the left side and that on the right side 
have opposite fermion number parity. When the per- 
pendicular vortex lines cross the two right-handed chi- 
ral vortex lines, a similar level crossing occurs for the 
right-handed fermions. Therefore the consequence of the 
instanton process is the pumping of a fermion from left- 
handed Fermi surface to the right-handed Fermi surface. 
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a b 




FIG. 2: (a) A system with two vertical vortex lines (lines 
with double arrow) and four horizontal chiral vortex lines 
(lines with single arrow). The two blue (red) horizontal lines 
are chiral vortices of Ol {9r), respectively (see text), (b) 
When one vertical vortex lines move across the two horizon- 
tal chiral vortex lines, two Majorana fermion zero modes are 
created along the horizontal vortex lines, (c) When both ver- 
tical vortex lines cross the chiral vortex lines, the Majorana 
fermions disappear, (d) The same process happens to the 
right-handed chiral vortex lines when the two vertical vortex 
lines are moved across them. Due to periodic boundary con- 
dition, the configuration in (a) and (d) are considered as the 
same. In term of energy spectrum, this process corresponds to 
the level crossing illustrated in (e), with a,b,c,d marking the 
time corresponding to the configurations showing in (a)-(d). 



B. Chiral vortex lines and s-wave proximity effect 

Since there is no independent U{1) symmetry for 8l 
and 6fi, one may expect the chiral vortex lines of 9l and 
6fj to be confined with each other, which seems to make it 
difficult to physically realize the vortex motion discussed 
above. However, in the following we will show that such 
chiral vortex lines can actually be realized experimentally 
on the interface between a topological superconductor 
and a trivial s-wave superconductor. By making use of 
the interface, the chiral vortex lines are deconfined with 
each other as long as they stay in the interface plane. 

As is discussed earlier, a trivial s-wave superconduc- 
tor with iV = corresponds to setting 6^ = 9^. If the 
phases are {9l, dn) = (0, 0) for the trivial superconductor 
and {9l,6b.) = (0, tt) for the nontrivial one, the interface 
between trivial and topological superconductors is a tt 
phase domain wall of dfj. (Another possibility is to 
have {OljSr) ~ (tt, 0) in the topological superconductor, 
so that the interface is a 9^ domain wall. This case can 
be analyzed in the same way.) If time- reversal symmetry 
is preserved, Or has to jump from to tt in the inter- 
face between the two materials, which means the pairing 
amplitude |Afl| must vanish somewhere near the surface. 
If time-reversal symmetry breaking is allowed on the in- 
terface, 9fi can change from to tt smoothly across the 



surface. Since 9^ is periodic in 27r, there are different 
ways to interpolate between = and 9fi = tt, corre- 
sponding to 

A6'fl = J dzdz9R = TT + 2mr 

A vortex of 9ii can be created if we consider a domain 
wall between different interpolations. If A6'fl = — tt for 
X < and A9b. = tt for a: > 0, a chiral vortex line is 
created at the line a; = 0, z = 0. 

Such a chiral vortex line can be realized in a tri- 
junction configuration shown in Fig. [3] which consists 
of two s-wave superconductor films on the surface of 
the 3D topological superconductor)^ The phases of the 
superconductors can be tuned such that (0L,^i?J = 
±(7r/2,7r/2) for the two s-wave superconductors, and 
i(^L,9ji) ~ (0,7r) for the topological superconductor. 
Across the junctions between the three superconductors, 
the phases 9l and fl/j interpolate continuously. There- 
fore if we follow the evolution of the phases along the 
counterclockwise circle around the junction (shown by 
the red circle in Fig. [3] (a)), the evolution of 9l and 9r 
are shown in Fig. [3] (b) by the blue and red paths, re- 
spectively. From this picture we see immediately that the 
winding number and Nn of 9]^, 9f{ are different and 
we always have ~ Nn — 1. Depending on whctlrer 
the interpolation from —tt/2 to +tt/2 across the junction 
between two s-wave superconductors is done through 
or TT, we may have {NljNh) = (0,-1) or (1,0), both 
of which correspond to the same Major ana- Weyl fermion 
propagating along the junction as is illustated in Fig. [3] 
(a) by the dashed line with arrow. 

The Majorana- Weyl fermion along the junction in this 
case can also be understood from an alternative surface 
state point of view. As has been discussed in Ref.— , 
the s-wave superconductor on the surface with phase 7r/2 
breaks time-reversal symmetry and leads to a mass term 
of the Majorana fermion on the surface. The surface 
state of the TSC is described by the low energy effective 
theory 

j (fxif [ax{-idx) + (yz{-idy) + m{x)<7y\ri (41) 

with •m{x) = m{—m) for a; > and x < 0, respectively. 
The Josephson junction is a mass domain wall for the 
surface Majorana fermion, which induces a Majorana- 
Weyl fermion bound state propagating on i^M^. 

This construction can be generalized to systems with 
multiple Fermi surfaces. For a TSC with 9i = Q or tt for 
each Fermi surface, the s-wave Josephson junction shown 
in Fig. [3] is a vortex line for each Fermi surface with 
9i = TT. If the Fermi surface has Chern number Cu, the 
chiral vortex line has Cu Majorana- Weyl fermion modes. 
(If Cu < 0, the corresponding Majorana- Weyl fermion 
has opposite chirality.) Therefore the total number of 
Majorana- Weyl fermions is 

N = -j:C'u'—f^ = lT.^ue^'' (42) 

i i 
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FIG. 3: (a) Illustration of the tri-junction configuration with 
two s-wave superconductors on top of a topological supercon- 
ductor (TSC). The phase of the s-wave superconductors are 
— 7r/2 in the blue region and 7r/2 in the yellow region. The 
vector (0,7r), (7r/2,7r/2) etc stands for the value of {9l,9r) 
in the corresponding regions, (b) The evolution of phases 8l 
(blue solid line) and 9r (red dashed line) along the red cir- 
cle path marked in panel (a). In such a configuration 9r has 
winding number —1 while 9l has winding number 0. Depend- 
ing on the sign of the phase gradient between the two s-wave 
superconductors (the position of the vertical path connect- 
ing 7r/2 and ~n/2 relative to the origin), one may also have 
winding number for 9r and 1 for Ol- Independent from this 
choice, there is always a Majorana-Weyl fermion propagating 
along the junction, as is illustrated by the dashed line in panel 
(a). 

which agrees with the bulk topological invariant of the 
TSC. 



V. CONCLUSION AND DISCUSSION 

In conclusion, in this work we have derived a new topo- 
logical field theory description of the three-dimensional 
TRI TSC. More precisely, our theory describes the dy- 
namics of electromagnetic field coupled with the dynamic 
phase fluctuation of the superconducting order parame- 
ter in a generic (3-|-l)-d superconductor. The topological 
property of the TSC is determined by the ground state 
value of the superconducting phases, so that the physical 
distinction between superconductors with different topo- 
logical invariants is characterized by our theory. The su- 
perconducting phases couple to electromagnetic field in 
an axionic topological coupling, with the coefficient de- 
termined by a topological property — the Chern number 
of each Fermi surface. As physical consequences of our 
topological field theory, we discussed the spectral flow 
and fermion number parity change in the process of mu- 
tual linking between chiral vortex lines. Physically, the 



chiral vortex lines can be realized by Josephson junctions 
between s-wave superconductor films on top of a 3d TRI 
TSC, so that the prediction of our theory can be verified 
in future experiments for possible candidate materials of 
3d TSC. 

Following the same derivation as that in Sec. IIIIl one 
can also investigate the response to gravitational field 
which was discussed in Refs i^^i'^^ . If we consider the cou- 
pling of fermions to gravitational field in the (4 + l)-d 
theory ([T^ . beside the Chern-Simons term one will 
obtain another mixed Chern-Simons term 

5grav = ^ I d'^e^''''A,Rl,^R%^ (43) 

with Riemann curvature. If we only consider the 

gravitational fiuctuation of the physical (3-|-l)-d manifold 
and keep the extra dimension fiat, R^^^^ is nonvanishcd 
only if no index is in the extra dimension. Therefore this 
topological term is reduced to 

^^^^^ = 1^3^ J d'^^'-'-'^RI^.RLr (44) 

This action can be straightforwardly generalized to the 
more generic form in analogy to Eq. ([24]): 

^^'■^^ = T5^E^i^^' / d'^e^^^^Rt^.RLr (45) 

i 

This result is consistent with the results of Ref<^i^. 
However, the apparent contradiction between the integer- 
valued topological invariant of the bulk TSC and the Z2 
value of the coefficient 9^ —Or = 0,7r in the presence 
of time-reversal symmetry is resolved in the new deriva- 
tion, since Ol^Oh are now interpreted as superconduct- 
ing phase fields rather than a topological order param- 
eter. The bulk topological invariant is determined by 
the integer coefficients Cu of this coupling together with 
the ground state value of the superconducting phase, as 
is given in Eq. (|25p. The main difference between the 
gravitational and electromagnetic response is that there 
is no Higgs term in the gravitational response, so that 
the gravitational field still propagates in the supercon- 
ductor (if there is an Einstein-Hilbert term in the ac- 
tion). Along a chiral vortex line oi di, a gravitational 
anomaly^ emerges which is consistent with the presence 
of a Majorana-Weyl fermion mode along it. 
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